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Abstract 

A technique for approximating the modal 
aerodynamic influence coefficients [AIC] 
matrices by using basis functions has been 
developed and validated. An application of the 
resulting approximated modal AIC matrix for a 
flutter analysis in transonic speed regime has 
been demonstrated. This methodology can be 
applied to the unsteady subsonic, transonic and 
supersonic aerodynamics. The method requires 
the unsteady aerodynamics in frequency- 
domain. The flutter solution can be found by 
the classic methods, such as rational function 
approximation, k, p-k, p, root-locus et cetera. 
The unsteady aeroelastic analysis for design 
optimization using unsteady transonic 
aerodynamic approximation is being 
demonstrated using the ZAERO™ flutter solver 
(ZONA Technology Incorporated, Scottsdale, 
Arizona). The technique presented has been 
shown to offer consistent flutter speed 
prediction on an aerostructures test wing 
[ATW] 2 configuration with negligible loss in 
precision in transonic speed regime. These 
results may have practical significance in the 
analysis of aircraft aeroelastic calculation and 
could lead to a more efficient design 
optimization cycle. 

1 Introduction 

Supporting the Aeronautics Research Mission 
Directorate guidelines, the National 
Aeronautics and Space Administration [NASA] 
Dryden Flight Research Center [DFRC] is 
developing an object-oriented multidisciplinary 
design, analysis, and optimization [MDAO] 


tool [1]. This tool will leverage existing tools 
and practices, and allow the easy integration 
and adoption of new state-of-the-art software. 
At the heart of the object-oriented MDAO tool 
is the central executive module [CEM] as 
shown in figure 1 . In this module, the user will 
choose an optimization methodology; provide 
starting and side constraints for continuous as 
well as discrete design variables and external 
file names for performance indices, which 
communicate between the CEM and each 
analysis module; submit script commands to 
prepare input data for each analysis code; 
execute analyses codes; compute performance 
indices using post-processor codes; and 
compute an objective function and constraint 
values from performance indices. The 
structural analyses modules such as 
computations of the structural weight, stress, 
deflection, buckling, and flutter and divergence 
speeds have been developed and incorporated 
into the object-oriented MDAO framework. 

With the increased complexity of the 
configuration and the need to increase the 
fidelity of the nonlinear aerodynamic equations 
to be solved in transonic regime, the 
computational cost of the unsteady aeroelastic 
analysis required to generate aerodynamic 
influence coefficients [AIC] matrices for full 
aeroelastic response and design optimizations 
is very high [1], Also, these analyses are 
usually performed repeatedly to optimize the 
final design. Even though the computational 
cost may be reduced by the use of advanced 
algorithms and improved computer hardware 
processing speeds, these full aeroelasticity 
analyses cannot be incorporated effectively 
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within a preliminary design and optimization 
environment. 

The integration of traditional unsteady 
computational fluid dynamics [CFD] into 
MDAO involving disciplines such as 
aeroelasticity and aeroservoelasticity is, at 
present, a costly and impractical venture [1]. At 
industry, aircraft designs using nonlinear CFD 
codes are mainly based on the manual trial and 
error approach, and this will considerably slow 
down the whole design procedure. As a result, 
there is considerable motivation to be able to 
perform unsteady aeroelastic calculations more 
quickly and inexpensively. The development of 
a MDAO tool to automate this process is 
needed to reduce time required for the early 
conceptual and/or preliminary design stages. 

The results of the rapid computation of a 
modal AIC matrix and aeroelastic response 
calculations using an approximated modal AIC 
matrix based on basis function approximation 
[BFA] has been previously validated and 
demonstrated successfully in subsonic speed 
regime to reduce the computation time. This 
approach reduces the modal AIC computation 
time using existing AIC matrices during 
optimization procedure by 80 percent [2], This 
paper extends the study and validation of the 
application of the BFA from the subsonic speed 
regime to the transonic speed regime. 

The transonic speed range introduces some 
particular flutter problems that are not 
commonly treated by classic linear flutter 
theory. Classic coupled torsion and bending 
wing flutter occurs when, with increase in 
dynamic pressure, the bending frequency 
increases and the torsional frequency decreases 
until the separation between these two 
frequencies becomes sufficiently small to allow 
exchange of energy between the two 
corresponding modes of the structure. In 
contrast, for wings with attached flow at 
transonic speeds, the flutter mechanism is 
somewhat different. It is very important to 
determine the location of the shocks, which is 
known to have fundamental importance on the 
correct estimation of the unsteady aerodynamic 
loads, for the transonic flutter prediction. The 
capture of the shock is impossible using any 


linear code at transonic speed. The general 
aeroelastic problem involving the interaction 
between the aircraft structure and the unsteady 
transonic flow with the presence of shock 
waves can only be solved using CFD codes, 
based on transonic small disturbance [TSD] 
theory, Euler codes, or the complete Navier- 
Stokes equations. The order of complexity and 
computational costs increases from the level of 
the TSD theory to the Navier-Stokes theory 
that can accurately predict interactions between 
boundary layers and local shocks on typical 
aerodynamic surfaces. In this paper, CAP- 
TSDv [3] code is used to capture the shock at 
steady state and flutter prediction. 

This paper presents the flutter results of 
the aerostructures test wing [ATW] 2 model 
without tip boom, in transonic speeds using 
BFA and traditional methods also called the 
direct method throughout the rest of the paper. 
The MSC NASTRAN (MSC Software 
Corporation, Santa Ana, California) [4], 
CAP-TSDv, ZAERO ZONA6 linear code, and 
ZTAIC nonlinear code [5] (ZONA Technology 
Incorporated, Scottsdale, Arizona) are used in 
this study for computing the flutter boundary. 

2 Basis Function Approximation 

In linear algebra, a basis is a set of vectors that, 
in a linear combination, can represent every 
vector in a given vector space or free module, 
and such that no element of the set can be 
represented as a linear combination of the 
others [6]. In other words, a basis is a linearly 
independent spanning set. Since mode shapes 
are orthogonal to each other and each mode has 
a natural frequency associated with the 
eigenvalue, structural mode shapes obtained 
from different mass or various stiffness 
configurations of the airplane can be used as 
the basis functions. These basis functions are 
comparison functions [7] since all the 
geometric and natural boundary conditions of 
the airplane are satisfied. The mode shapes of 
the airplane with any of the configurations in 
the design spaces can be used to generate the 
basis functions. Furthermore, each mode shape 
of the target configuration of the airplane can 
be approximated as a linear combination of a 
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set of basis functions. This basis function 
approximation has been described in reference 
2 . 

In general, a time-domain approach is 
more convenient for analysis tasks and 
nonlinear problems; on the other hand, a 
frequency-domain approach is better for design 
tasks and linear problems. Most modern CFD 
codes were developed using the time-domain 
approach, and therefore, the results of these 
CFD runs are usually the time-histories of 
pressure distributions, structural displacements 
et cetera. To use the modem CFD computation 
for our design optimization studies, we are 
interested in developing an approximate 
unsteady computational aeroelastic analysis in 
the frequency-domain to dramatically reduce 
the required computation time for optimization 
runs in the transonic speed regime. To 
accomplish these objectives, the BFA method 
in reference 2 can be extended and evaluated 
for approximation in transonic speeds. The 
BFA method requires the unsteady 
aerodynamics to be represented in the 
frequency domain. The flutter solution can be 
found by classic methods, such as rational 
function approximation, k, p-k, p, g method, 
and root-locus et cetera. 

The computation of a modal AIC matrix 
and aeroelastic response calculations in 
transonic speeds will be discussed in this paper 
and the approximation process is outlined in 
the flowchart given in figure 2. 

In step 1, a set of representative basis 
functions, T , is chosen with the intent to 
capture salient features of the modal responses 
the airplane is expected to encounter in the 
various design spaces. 

In step 2, the steady state pressure 
distributions for the transonic speed analysis 
can be obtained by using a CFD code or wind 
tunnel test. 

In step 3, a set of representative basis 
modal AIC matrices, Q , are computed 
corresponding to the representative basis 
functions defined in step 1 and steady pressure 
distributions obtained in step 2 at any Mach 
number and reduced frequencies. These basis 
modal AIC matrices are used as input for 


approximate modal AIC matrix calculation in 
step 5. These basis modal AIC matrices are 
computed only once and stored before starting 
optimization. 

In step 4, during optimization, for a set of 
given structural mode shapes, 0 , each mode 
shape is decomposed in a linear combination of 
the basis functions, V F . The i-th mode shape, 
0,, (i=l, 2, ..., m, where m is the number of 
mode shapes) is approximated through the use 
of a least squares fitting together with the 
following basis functions as shown in 
equation (1), 

h-tti'i'k (i) 

k = 1 

where \f/ k is the k-th basis function and a 
coefficient, [3 ’ k , is the modal participation 
factor of the k-th basis function on the i-th 
mode shape and n is the number of basis 
functions. 

In step 5, an approximate modal AIC 
matrix, Q , is computed based on a basis modal 
AIC matrix, Q , and modal participation factors 
in equation (2), 

Q,i=Qu=iim&r ( 2 > 

r-\s=\ 

where Q is the i-th row and j-th column 

element of the modal AIC matrix Q. 

In step 6, the flutter analysis is performed 
using the approximate modal AIC matrix in 
equation (2). The flutter speed and frequency 
could be used for optimization. 

3 Applications 

In an effort to validate the BFA in transonic 
speed, a modified ATW2 model was chosen 
and the flutter results using approximate 
method were compared with the direct flutter 
results for both subsonic and transonic speed. 
To simplify the problem, especially for CFD 
mesh generation in CAP-TSDv, the wing tip 
boom is removed from this application. The 
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original ATW2 [8] with wing tip boom is the 
test article actually designed, built, and tested 
at NASA DFRC as shown in figure 3. The 
ATW2 was clamped onto a circular plate, 
which was bolted to a mounting panel. The 
mounting panel was attached to the flight test 
fixture [FTF] [9]. The FTF was installed on the 
centerline F-15B (McDonnell Douglas, now 
The Boeing Company, Chicago, Illinois) pylon 
in flight configuration. 

For the direct method, the flutter analysis 
of the target design configuration is performed 
using the ZAERO linear subsonic method 
(ZONA6) and nonlinear transonic method 
(ZTAIC) with the first 10 structural modes and 
16 reduced frequencies at 4 Mach numbers 
(Mach 0.80, 0.90, 0.95, and 0.98). The steady 
state pressure input data as the initial condition 
for the ZTAIC method used in the present 
analysis are obtained by using the CAP-TSDv 
code at all four Mach numbers. In order to 
prove the setup and the process of the ZTAIC 
method correctly, the flutter results from both 
the ZONA6 and ZTAIC methods are compared 
and used to verify the process at subsonic 
Mach numbers. At a subsonic Mach number 
like Mach 0.80, before any shock is formed, 
both the linear and nonlinear solver should 
predict similar flutter results. Once the process 
was verified, the ZTAIC method was used for 
the flutter prediction in the transonic speed 
regime. 

The CAP-TSDv calculations were 
performed on a 150x80x80 points 
computational grid with 41 points along each 
of the 57 spanwise chords on the wing as 
shown in figure 4. The steady pressure 
distributions obtained from the CAP-TSDv 
computations at Mach 0.80, 0.90, 0.95, and 
0.98 all at 0.0 deg and 1.0 deg of angle of 
attack [AO A] are shown in figure 5 and 6; and 
the steady pressure differences between the 
upper and the lower surfaces under 1.0 deg 
AOA at the same Mach numbers are shown in 
figure 7. With 3 percent structural damping, the 
primary flutter speed and flutter frequency 
comparison of the ZONA6 and ZTAIC 
methods are summarized in figure 8. At the 
subsonic Mach number of 0.80, the ZTAIC 


results are in close agreement with those of the 
ZONA6, as expected, since a transonic shock 
has not formed yet. 

The primary flutter mode is the first 
bending mode at Mach 0.80. However, the 
second bending mode becomes the primary 
flutter mode when the Mach number is 0.98. 
On the other hand, the flutter mode interchange 
does not happen with the ZONA6 computation. 
At Mach 0.98, the primary flutter mode with 
the ZTAIC computation is the second mode, 
81 Hz, as shown in figure 8, while the primary 
flutter mode with the ZONA6 computation is 
the first mode, 55 Hz. This mode interchange 
may be caused by the local high pressure zone 
near the center of the wing as shown in figures 
5 and 6. 

Based on the pressure differences in figure 
7, the aerodynamic load distributions at 
Mach 0.80 and 0.90 are fairly similar in shape, 
which are basically the subsonic distributions, 
(i.e. high near the leading-edge area and low 
near the trailing-edge area). When the Mach 
number became 0.95, the area of the high 
aerodynamic loads was extended to the aft 
section of the wing-root chord, and an even 
higher load area was formed near the leading- 
edge of the wing mid chord. At Mach 0.98, 
aerodynamic load distributions are changed 
drastically, and compared to the Mach 0.95 
case, highly concentrated aerodynamic load 
distributions are observed from the leading 
edge of the wing tip chord to the trailing edge 
of the wing root chord. From these 
observations, we may conclude that the 
transonic effects start to happen between 
Mach numbers of 0.90 and 0.95, and strong 
shocks are formed at Mach 0.98 as shown in 
figures 5 and 7. Probably the second bending 
mode in figure 9 can be easily excited by the 
concentrated aerodynamic load distributions at 
Mach 0.98. 

In this ATW2 example, five different 
design configuration points in the parameter 
design space are selected for generating basis 
functions for the approximate flutter solution. 
The design variations used in this study are 
variations in wing thickness and stiffness that 
are effected by variations in the number and 
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orientation of the composite laminate plies. 
These design variations are intended to capture 
salient features of the modal responses 
resulting from mass and stiffness variations in 
the design space. The structural mode shapes 
obtained from those five design configurations 
are used as the basis functions. The finite 
element model and the first three mode shapes 
of the target configuration are shown in figure 
10. The total weight and center of gravity of 
the various fictitious and target design 
configurations are listed in table 1. The target 
configuration in this paper is simply an 
arbitrary design point within the design space 
and is used as a check case. 

First, 10 structural mode shapes for each 
design configuration (a total of 50 basis 
functions) are generated using 
MSC/NASTRAN. The natural frequencies of 
the various fictitious and target design 
configurations are listed in table 2. Next, the 
size of 50 by 50 basis modal AIC matrix 
corresponding to those 50 basis functions is 
computed for each of the 16 reduced 
frequencies at Mach 0.80, 0.90, 0.95, and 0.98 
using ZONA6 and ZTAIC. To solve the modal 
participation factors, the mode shapes of the 
target design configuration are fitted using the 
basis functions with a least squares method. 
Then the approximate modal AIC matrices are 
computed based on the modal participation 
factors. Finally, the approximate flutter 
boundaries can be solved by using any linear 
frequency domain flutter solution code and the 
approximated modal AIC matrices. It is 
important to point out that in order to provide 
meaningful comparisons of the direct method 
with the BFA method, the ZTAIC method is 
used to predict the flutter boundaries in this 
paper. Given that the BFA approach can be 
used with any kind of CFD code, ZTAIC is 
selected for CFD simulation in this paper due 
to limited resources. The focus in this paper is 
to demonstrate that the BFA method can be 
used to predict adequate flutter speed and 
flutter frequency at transonic speed regime 
regardless the level of CFD code. 

As shown in table 3, using ZTAIC, both 
methods predicted similar flutter speed and 
flutter frequency at Mach 0.80, 0.90, 0.95, and 


0.98. The speed versus damping, V-g, and 
speed versus frequency, V-co, curves from the 
direct and BFA methods at Mach 0.90, 0.95, 
and 0.98 are given in figures 11, 12, and 13. 
These figures indicate that an interchange 
between the first bending and the second 
bending flutter modes results when Mach 
number increased from 0.90 to 0.98. At 
Mach 0.90, the first bending mode, 56 Hz, is 
predicted as the primary flutter mode. This 
flutter mode became the hump mode at 
Mach 0.95 and stable mode at Mach 0.98. 
Meanwhile, the stable second mode at 
Mach 0.90 became unstable at Mach 0.95 and 
later became primary flutter mode at 
Mach 0.98. With 3 percent structural damping, 
the primary flutter speed and flutter frequency 
are summarized in table 3 and figure 14. The 
largest percentage error of flutter speed and 
frequency is 1.60 percent at Mach 0.98 and 
1.27 percent at Mach 0.80 respectively. 
Although a different CFD code may generate 
different unsteady AIC and predict different 
flutter results, the important point is that the 
results of the BFA method is following the one 
predicted by the direct method. Overall, the 
resulting flutter boundaries prediction using the 
BFA method is well compared with the direct 
method in transonic speed. The flutter 
boundaries prediction based on the BFA 
method are essentially identical to those using 
direct solution on the target design 
configuration. 

Table 4 lists the computational cost of 
computing the modal AIC matrices with given 
structural mode shapes at one Mach number 
and 16 reduced frequencies using the direct 
method and the BFA method. The size of the 
ATW2 model is about 260 degrees of freedom. 
The computational cost of generating the 
modal AIC using the direct method and the 
BFA method is about 43 s and 8 s using 
ZONA6 and 49 s and 8 s using ZTAIC 
respectively in this ATW2 example. As a 
result, the BFA approach reduces the modal 
AIC computation time using existing AIC 
matrices by 80 percent. The comparison was 
done on an Intel® Core™2 2.80GHz CPU 
computer (Intel Corporation, Santa Clara, 
California). The computational cost of 
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generating the basis functions before the 
parametric study is not included in the 
comparison. 

4 Concluding Remarks 

Many of the most interesting engineering 
applications are very complex systems, such as 
aeroelastic analyses and optimization in which 
the analysis is achieved through the use of 
expensive and time consuming numerical 
simulations. The processing times on even the 
most powerful workstations may be too 
expensive. A technique for approximating the 
modal AIC matrix by using basis function 
approximation has been applied in transonic 
speed regime, and a process for using the 
resulting AIC matrix in aeroelastic analyses has 
been developed and demonstrated. The 
technique presented has been shown to offer 
consistent flutter speed prediction on an ATW2 
configuration with insignificant small loss in 
accuracy, less than 1.6 percent error in the 
ATW2 test case, and the results are essentially 
identical to those using direct solution. The 
BFA approach has been demonstrated 
successfully in the transonic speed regime to 
reduce the computation time of generating 
modal AIC using basis functions. The BFA 
approach has shown an 80 percent reduction in 
the modal AIC computation time using existing 
AIC matrices. These results may have practical 
significance in the aeroelastic analysis of an 
aircraft and could lead to a more efficient 
design optimization cycle in a transonic speed 
regime. For the design optimization and 
parameter study, hundreds or even thousands of 
analyses have to be performed. Creation of an 


approximated modal AIC matrix using basis 
functions allows for analysis, and a parameter 
optimization cycle of the system in minutes, 
rather than hours of computer time on common 
personal computers. 
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Tables 

Table 1. Summary of total weight and center of gravity location for the ATW2 without the tip boom model with different 
fictitious design configurations. 


Fictitious configuration 

Config. 1 

Config. 2 

Target 

Config. 3 

Config. 4 

Config. 5 

Total weight, lb 

2.38 

2.53 

2.81 

2.88 

3.66 

4.09 

X C g, in 

12.56 

12.70 

12.93 

12.98 

13.41 

13.58 

Y C g, in 

-7.17 

-7.24 

-7.35 

-7.38 

-7.59 

-7.67 

Z C g, in 

-3.12 

-2.95 

-2.66 

-2.59 

-2.05 

-1.85 


Table 2. Summary of natural frequencies (Hz) of the ATW2 without the tip boom model with different design 
configurations. 


Mode 

Config. 1 

Config. 2 

Target 

Config. 3 

Config. 4 

Config. 5 

1 

23.08 

22.77 

22.06 

21.88 

20.01 

19.13 

2 

85.61 

85.76 

84.85 

84.49 

79.34 

76.50 

3 

120.67 

126.71 

131.30 

131.62 

126.99 

122.32 

4 

147.99 

168.97 

194.28 

194.17 

183.32 

177.44 

5 

193.02 

195.30 

194.80 

198.63 

223.60 

228.76 

6 

261.81 

277.74 

291.27 

292.60 

286.04 

277.39 

7 

315.54 

322.42 

323.05 

322.07 

304.26 

295.76 

8 

399.82 

423.26 

442.11 

443.74 

434.38 

426.50 

9 

459.62 

472.40 

472.64 

470.78 

442.33 

431.65 

10 

560.99 

582.69 

595.67 

595.64 

572.74 

561.21 


Table 3. Summary of the ATW2 without the tip boom model flutter results comparison for BFA and the direct method 
using ZTAIC. 


Mach 

Flutter speed, KEAS 

Flutter frequency, Hz 

BFA 

Direct 

Error 

BFA 

Direct 

Error 

0.80 

621.77 

622.94 

-0.19% 

59.85 

60.62 

-1.27% 

0.90 

612.39 

612.38 

0.00% 

56.56 

56.56 

0.00% 

0.95 

641.94 

640.68 

0.20% 

53.63 

54.01 

-0.70% 

0.98 

557.56 

566.60 

-1.60% 

81.23 

81.10 

0.16% 


Table 4. Computational cost comparison of BFA and the direct method. 



ZONA6 

ZTAIC 

Direct method elapsed time 

43 s 

49 s 

BFA method elapsed time 

8 s 

8s 
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Figures 
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Fig. 1. Object-oriented multidisciplinary design, Fig. 4. CAP-TSD wing surface mesh, 

analysis, and optimization tool. 



Fig. 2. Flowchart of the flutter analysis module in the 
object-oriented MDAO tool. 
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Fig. 5. Steady pressure distributions on the upper surface 
using CAP-TSDv code (0 deg angle of attack). 
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Fig. 6. Steady pressure distributions on the upper surface 
using CAP-TSDv code (1 deg angle of attack). 


Fig. 3. Aerostructures test wing 2. 
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Fig. 7. Steady pressure difference under 1 deg angle of 
attack using CAP-TSDv code. 
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Fig. 8a. Flutter speed. 
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Fig. 8b. Flutter frequency. 

Fig. 8. Flutter results for the ATW2 without tip boom 
(Z0NA6 compared with ZTAIC). 
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Fig. 9. V-g and V-co plots for the ATW2 at Mach 0.98 
using Z0NA6 and ZTAIC methods. 
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Fig. 10. The aero structures test wing 2 without tip boom 
mode shapes. 



no no «o »o «» a» m m 300 ooo soo ao no mo 

Speed (Meat) Speed (Keai) 


Fig. 11. V-g and V-co plots for the ATW2 at Mach 0.90 
using direct ZTAIC and BFA ZTAIC methods. 
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Fig. 12. V-g and V-co plots for the ATW2 at Mach 0.95 
using direct ZTAIC and BFA ZTAIC methods. 
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Fig. 14b. Flutter frequency. 

Fig. 14. Flutter results for the ATW2 without tip boom. 
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Fig. 13. V-g and V-co plots for the ATW2 at Mach 0.98 
using direct ZTAIC and BFA ZTAIC methods. 
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Introductions 


□ In support of the Aeronautics Research Mission Directorate (ARMD) guidelines, the 
NASA Dry den has developed an Object-Oriented Optimization (O 3 ) tool. This tool 
leverages existing tools and practices, and allow the easy integration and adoption 
of new state-of-the-art software. 


❖ 


❖ 


❖ 


❖ 


❖ 


Choose an optimization method. 

Define objective and constraint functions 
from performance indices. 

Provide starting and side constraints for 
continuous as well as discrete design 
variables and executable and I/O file 
names for each module. 

Submit script commands for executing 
each module. 

Structural analyses modules such as 
computations of the structural weight, 
stress, deflection, buckling, and flutter 
and divergence speeds as well as a 
structural model update module have 
been developed and incorporated. 


Need to develop; 
Depends on 
Analysis code 



Object-GiientedMDAG tool 
with model tuning capability' 


Wesley Li-3 



Introductions (continued) 


□ Unsteady aeroelastic analysis and optimization is a challenging task 


❖ With the increased complexity of the model configuration and the need to increase 
the fidelity of the nonlinear aerodynamic equations, huge DOF problems, to be 
solved, the computational cost of the unsteady aeroelastic analysis required to 
generate aerodynamic influence coefficients (AIC) matrices for full aeroelastic 
response is very high, especially in transonic regime. 

❖ Also, these analyses are usually performed repeatedly to optimize the final design. 
Even though the computational cost may be reduced by the use of advanced 
algorithms and improved computer hardware processing speeds, these full 
aeroelasticity analyses cannot be incorporated effective within a preliminary 
design and optimization environment. 

❖ The integration of traditional unsteady computational fluid dynamics (CFD) into 
MDAO involving disciplines such as aeroelasticity and aeroservoelasticity is, at 
present, a costly and impractical venture. At industry, aircraft designs using 
nonlinear CFD codes are mainly based on the manual trial and error approach, 
and this will considerably slow down the whole design procedure. 
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Objectives 


□ To be able to perform aeroelastic calculations more quickly 
and inexpensively by reducing the computation time for 
generating modal AIC matrices using Basis Function 
Approximation (BFA) method during the optimization 
procedure. 
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Basis Function 


□ In linear algebra, a basis is a set of vectors that, in a linear combination, can represent every vector 
in a given vector space or free module, and such that no element of the set can be represented as a 
linear combination of the others. 


□ Mode shapes are orthogonal to each other and each mode has a natural frequency associated with 
the eigenvalue, therefore, structural mode shapes can be used as the basis functions. 

□ Mode shapes are comparison functions for the original eigenvalue problem since all the geometric 
and natural boundary conditions of the airplane are satisfied. 


□ Mode shapes obtains from different mass or various 
stiffness configurations of the airplane can be used to 
capture significant features of the modal responses of the 
airplane to encounter in the various design space. 


Mode shapes are Comparison Functions 


CQ 
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Natural B.C. 




Basis Function Approximation (BFA) 


□ The results of the rapid computation of a modal AIC matrix and aeroelastic response calculations 
using an approximated modal AIC matrix based on basis function approximation (BFA) has been 
previously validated and demonstrated successfully in subsonic speed regime to reduce the 
computation time. (AIAA-201 0-3085, pp 1-18, April 12-15, 2010) 

Compute unsteady aerodynamics 
in the frequency domain using 
Basis Function Approximation 
method. 

Approximate Modal AIC has been 
validated by comparing the 
aeroelastic solution instead of just 
comparing each of the elements in 
the matrices. The flutter solution 
can be found by the classic 
methods, such as rational function 
approximation, k, p-k, p, g- 
method etc. 

The BFA approach reduces the 
modal AIC computation time 
using existing AIC matrices 
during optimization procedure by 

80 percent. 


Direct During Optimization: Unchanged During Optimization: Changed 

Modal AIC — I ndof: number of Degrees of Freedom 

^ij 



General AIC: ndof x ndof 


Read matrix A: expensive 


Direct vs. Basis Function Approximation 


BFA 


n 



n n 


Approximate / 
Modal AIC \J 


V 



ndof x n 

Least Squares Method 

Read matrices ip & Q : cheaper 
n x n 

n: number of Basis Function 


r = 1 s=\ 


□ This paper extends the study and validation of the application of the BFA from the 
subsonic speed regime to the transonic speed regime. 
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Flutter Analysis Module Flowchart 


□ In step 1 , generate a set of 
“representative” structural basis 
function using structural mode 
shapes. 

□ In step 2, compute the steady state 
pressure for CFD run 

□ In step 3, compute the basis modal 
AIC matrices corresponding to the 
structural basis functions defined 
in step 1 at each Mach number 
and reduced frequencies. These 
basis modal AIC matrices are used 
as input for approximate modal 
AIC matrix calculation in step 5. 


Database: generated before optimization 


Trs 


Generate basis 
functions 


Given structural 
mode shape 



4, 

-► 


i 


2 

Obtains steady 

ionic ? 

pressure 

Y 

(Aero) 

■ 

N 

\ 

3 

Generate basis 


AIC (Aero) 


T 


Fit mode shapes 
using basis 
functions 


A 


Compute 
approximate 
modal AIC 


6 

Flutter 


anal 

lysis 




Flutter speed 
& frequency 



Q 


Flutter analysis module 


Flutter Analysis Module Flowchart 
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Flutter Analysis Module Flowchart (continued) 


□ In step 4, decompose each given structural 
mode shape in a linear combination of the basis 
functions W using least squares fitting. The i-th 
mode shape, & i9 is defined as 


□ 


<t>, 


Pitt’ 


In step 5, compute the 
approximate modal AIC matrix 
based on a basis modal AIC 
matrices and modal 
participation factors. This 
approximate modal AIC matrix 
will be used for flutter analysis. 


n n 


Q» = 



tW,Q, 


r = 1 s= 


Given stinctural 
mode sli up e 


Database: generated before optimization 
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, , , 7 , Flutter Analysis Module Flowchart 

where n is the number of basis junctions, 

where Q tJ is the i-th row and j-th column element of the modal AIC matrix Q. 

where t// k is the k-th basis function and a coefficient, , is the modal participation factor of the k th basis 

function on the i-th mode shape and n is the number of basis functions. 
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Validations 



NASA 


□ 

□ 

□ 

□ 


In order to validate the BFA method in transonic flight regimes, the results of the flutter solution 
using this approximate method are compared with the direct flutter solution for both subsonic and 
transonic speed (analytical only). 

Test case: ATW2 without tip boom model. 


The original ATW2 with wing 
tip boom is the test article 
actually designed, built, and 
tested at NASA DFRC. 

TheATW2 was clamped onto a 
circular plate, which was 
bolted to a mounting panel. The 
mounting panel was attached to 
the flight test fixture (FTF). The 
FTF was installed on the 
centerline F-15B pylon in flight 
configuration. 


Centerline pylon 




Flight Test Fixture 


Mounting Panel 






ATW2 




Circular Plate 
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Validations (continue) 


□ To demonstrate that the BFA method can be used to predict adequate flutter speed and flutter 
frequency at transonic speed regime regardless the level of CFD code. 


□ To provide meaningful comparisons of the direct method with the BFA method, the identical flutter 
solver is used to predict the flutter boundaries. 


□ The BFA approach can be used with any kind of CFD code, ZAERO nonlinear transonic method, 
ZTAIC is selected for CFD simulation in this study due to limited resources. 


□ To prove the setup and the process of using ZTAIC method correctly, the flutter results from both 
the linear subsonic method, ZONA6 and ZTAIC methods are compared at subsonic Mach numbers. 

□ The flutter analysis of the target design configuration is performed with the first 10 structural modes 
and 16 reduced frequencies at 4 Mach numbers (Mach 0.80, 0.90, 0.95, and 0.98). 

□ The steady state pressure input data as the initial condition for the ZTAIC method used in the flutter 
analysis are obtained by using the CAP-TSDv code. 
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CAP-TSDv Steady Pressure Results 


□ The CAP-TSDv calculations were performed on 
a 150 x 80 x 80 points computational grid with 
41 points along each of the 57 spanwise chords 





Steady pressure difference under 1 deg angle of 
attack using CAP-TSDv code. 
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Steady pressure distributions on the upper surface 


using CAP-TSDv code (0 deg angle of attack). 
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Steady pressure distributions on the upper surface 
using CAP-TSDv code (1 deg angle of attack). 
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Flutter Results - Z0NA6 vs. ZTAIC 


□ At the subsonic Mach number of 0.80, the ZTAIC 
results are in close agreement with those of the 
ZONA6, as expected, since a transonic shock has not 
formed yet. 

□ At the transonic Mach number of 0.95 and 0.98, the 
flutter mode interchange does not happen with the 
ZONA6 computation. This mode interchange may be 
caused by the local high pressure zone near the center 
of the wing 
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A modified ATW2 Fictitious Configuration 

□ To generating basis functions for the approximation modal AIC and the approximate flutter 
solution, five different design configuration points in the parameter design space are selected. 

□ The design variations used in this study are variations in wing thickness and stiffness that are 
effected by variations in the number and orientation of the composite laminate plies. 

□ The design variations are intended to capture salient features of the modal responses resulting 
from mass and stiffness variations in the design space. 


Fictitious configuration 

Config. 1 

Config. 2 

Target 

Config. 3 

Config. 4 

Config. 5 

Total weight, lb 

2.38 

2.53 

2.81 

2.88 

3.66 

4.09 

X C g, in 

12.56 

12.70 

12.93 

12.98 

13.41 

13.58 

Y cg , in 

-7.17 

-7.24 

-7.35 

-7.38 

- 7.59 

-7.67 

^CG’ in 

-3.12 

-2.95 

-2.66 

-2.59 

-2.05 

-1.85 


Summary of total weight and center of gravity location for the ATW2 without the tip 
boom model with different fictitious design configurations (MSC NASTRAN). 
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Natural Frequencies & Mode Shapes 

□ Summary of natural 

frequencies of the ATW2 with 
different fictitious mass and 
stiffness configurations using 
MSC NASTRAN. 


□ Structural finite element model mode 
shapes of the ATW2 with target mass 
and stiffness configuration. 


Number of DOFs ~ 260 



Summary of natural frequencies (Hz) of the ATW2 without the tip 
boom model with different design configurations. 


Mode 

Config. 1 

Config. 2 

Target 

Config. 3 

Config. 4 

Config. 5 

1 

23.08 

22.77 

22.06 

21.88 

20.01 

19.13 

2 

85.61 

85.76 

84.85 

84.49 

79.34 

76.50 

3 

120.67 

126.71 

131.30 

131.62 

126.99 

122.32 

4 

147.99 

168.97 

194.28 

194.17 

183.32 

177.44 

5 

193.02 

195.30 

194.80 

198.63 

223.60 

228.76 

6 

261.81 

277.74 

291.27 

292.60 

286.04 

277.39 

7 

315.54 

322.42 

323.05 

322.07 

304.26 

295.76 

8 

399.82 

423.26 

442.11 

443.74 

434.38 

426.50 

9 

459.62 

472.40 

472.64 

470.78 

442.33 

431.65 

10 

560.99 

582.69 

595.67 

595.64 

572.74 

561.21 
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Flutter Results - Direct vs. BFA 


□ Summary of the ATW2 without the tip boom model flutter results comparison for BFA and the 
direct method using ZTAIC. 


Mach 

Flutter speed, KEAS 

Flutter frequency, Hz 

BFA 

Direct 

Error 

BFA 

Direct 

Error 

0.80 

621.77 

622.94 

-0.19% 

59.85 

60.62 

-1.27% 

0.90 

612.39 

612.38 

0.00% 

56.56 

56.56 

0.00% 

0.95 

641.94 

640.68 

0.20% 

53.63 

54.01 

-0.70% 

0.98 

557.56 

566.60 

-1.60% 

81.23 

81.10 

0.16% 
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Modified ATW2 Flutter Results - Vg Plots, Mach 0.90 


0.10 


V-g and V-co plots for the ATW2 at Mach 0.90 using direct ZTAIC and BFA ZTAIC methods. 
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Modified ATW2 Flutter Results - Vg Plots, Mach 0.95 



V-g and V-co plots for the ATW2 at Mach 0.95 using direct ZTAIC and BFA ZTAIC methods. 
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Modified ATW2 Flutter Results - Vg Plots, Mach 0.98 
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V-g and V-co plots for the ATW2 at Mach 0.98 using direct ZTAIC and BFA ZTAIC methods. 
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Computational Cost Comparison 


□ Computational cost comparison of BFA and Direct method. 

♦♦♦ The computational cost of computing the modal AIC with a given structural mode shape 
at one Mach number using the direct method and the BFA method. The comparison was 
done on an Intel® Core™2 2.80 GHz CPU computer. 

♦♦♦ The computational cost of generating the basis functions before the parametric study is 
not included in the comparison. 



ZONA6 

ZTAIC 

Direct method elapsed time 

43 sec 

49 sec 

BFA method elapsed time 

8 sec 

8 sec 


BFA is 5 times faster than Direct. 


□ For the design optimization and parameter study, hundreds or even thousands of analyses have 
to be performed. Creation of an approximated modal AIC matrix using basis functions allows 
for analysis, and a parameter optimization cycle of the system in minutes, rather than hours of 
computer time on common personal computers. 
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Concluding Remarks 


□ Many of the most interesting engineering applications are very complex systems, such as aeroelastic 
analyses and optimization in which the analysis is achieved through the use of expensive and time 
consuming numerical simulations. The processing times on even the most powerful workstations may 
be too expensive. 

□ The basis function approximation (BFA) has been validated and demonstrated successfully in the 
subsonic and transonic speed regime to reduce the computation time of generating modal AIC using 
basis functions. A process for using the resulting AIC matrix in aeroelastic analyses has been 
developed and demonstrated as well. 

□ The BFA technique presented has been shown to offer consistent flutter speed prediction on an ATW2 
configuration with insignificant small loss in accuracy, less than 1.6 percent error in the ATW2 test 
case, and the results are essentially identical to those using direct solution. The BFA approach has 
shown an 80 percent reduction in the modal AIC computation time using existing AIC matrices, 5 

times faster than direct method. 

□ For the design optimization and parameter study, hundreds or even thousands of analyses have to be 
performed. Creation of an approximated modal AIC matrix using basis functions may have practical 
significance in the aeroelastic analyses of an aircraft and could lead to a more efficient design 
optimization cycle. 
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Questions? 


